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Abstract 

The original Jaynes-Cummings model is described by a Hamiltonian which 
is hermitian and exactly solvable. Here we extend this model by several types 
of interactions leading to a non hermitian operator which doesn't satisfy the 
physical condition of space-time reflection symmetry (PT symmetry). The new 
Hamiltonians are either exactly solvable admitting an entirely real spectrum or 
quasi exactly solvable with a real algebraic part of their spectrum. 
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1 Introduction 



Several new theoretical aspects of quantum mechanics have been developped in the 
last years. In series of papers (see e.g. [1, 2] and [3] for a recent review) it is shown 
that the traditional self adjointness requirement of the Hamiltonian operator is not a 
necessary condition to guarantee a real spectrum and that the weaker condition of PT- 
invariance of the Hamiltonian is sufficient for that purpose. An alternative possibility 
for an operator to admit a real spectrum is also developed in [4]. It is the notion of 
pseudo- hermit icity. Following the ideas of [4], we remind here that a Hamiltonian is 
called T) pseudo-hermitian if it satisfies the relation r\Hrf x = H\ where r\ denotes 
a linear hermitian operator. It is this new notion (i.e pseudo-hermiticity property) 
of non hermitian Hamiltonians which explains the reality of their energy spectrum. 
This important property has further been considered in Refs. [5, 6]. 

Another direction of development of quantum mechanics is the notion of quasi 
exact solvability [7, 8] . It provides techniques to construct linear operators preserving 
a finite dimensional subspace V of the Hilbert space. Accordingly, the so called Quasi 
Exactly Solvable operators, once restricted on V can be diagonalized by means of 
algebraic methods. The QES property is strongly connected to finite dimensional 
representation of Lie or graded Lie algebras [7, 9, 10]. Amongst many models used to 
describe quantum properties of physical systems, the Jaynes-Cummings model play 
an important role [11, 12, 13, 14]. It describes, in a simple way the interaction of 
photons with a spin-1/2 particle. From the mathematical point of view, the Jaynes- 
Cummings model is described by a self-adjoint operator and it is completely solvable 
in a sense that the entire spectrum can be computed algebraically. 

The purpose of this paper is to consider operators generalizing the Jaynes-Cummings 
Hamiltonians which are neither self- adjoint nor PT-invariant but which are pseudo- 
hermitian with respect to two different operators. In particular, from the original 
Jaynes-Cummings model (JCM in the following), we construct an extended one by 
adding a polynomial of the form P{o),a) (a\a are the usual creation and annihi- 
lation operators) of degree d > 2 in the diagonal part of the hamiltonian. Some 
particular choices of P are constructed in such a way that the resulting operator be- 
comes QES. The non-diagonal interaction part is also modified in such a way that (i) 
multiple photon exchanges are allowed and (ii) the full operator can be hermitian or 
pseudo-hermitian. 

Here is the plan of the paper. In section 2, we revisit the Hamiltonian considered 
in Ref. [5] and express it in terms of differential operator of a real variable x. This 
reveals its exact solvability in terms of differential operators acting on sets of polyno- 
mials of appropriate degrees in x. In Sect. 3 we propose a family of operators which 
generalize the original JC Hamiltonian in several respects. The (pseudo)-hermiticity 
of these operators are analysed and the spectra and the eigenvectors are computed in 
details for a few of them. The differences in the spectrum corresponding to Hermi- 
tian and pseudo- Hermitian are pointed out. In particular, the energy eigenvalues are 
entirely real in spite of the fact that they are associated to a non hermitian and non 
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PT-invariant Hamiltonian. The reality of those eigenvalues is ensured by the pseudo- 
hermiticity of the Hamiltonians. The section 4 is devoted to QES extensions of the 
JCM. These are constructed in such a way that, both, one-photon and two-photons 
exchange terms coexist in the non-diagonal interacting terms. By construction, these 
new models preserve finite dimensional vector spaces of the Hilbert spaces ,the alge- 
braic part of the spectrum is computed in Sect 5. Further properties of these new 
types of QES operators, say H T , can be discussed. Namely, following the ideas of [15] 
we show in Sect. 6 that the solutions of the spectral equation H T tp = Et/j for generic 
values of E lead to new types of recurence relations. The relations between H T and 
specific graded algebras are pointed out in Sect 7. Finally, the section 8 is kept for 
concluding remarks. 



2 Exactly solvable pseudo-hermitian Hamiltonian 

In this section we consider the Hamiltonian describing a system of a spin-| particle in 

— * 

the external magnetic field, B which couples to a harmonic oscillator through some 
nonhermitian interaction [5] 

H = n& ■ B + hua^a + p(a + a — <7_a'). (1) 

Here a denotes Pauli matrices, p is some arbitrary real parameter and a± = \[a x ±.ia y }. 
a + and <7_ can be expressed in matrix form 



Our purpose is to relate the Hamiltonian above to an appropriate differential operator 
preserving a family of spaces of polynomials in the variable x, following the ideas of 
exactly and quasi-exactly solvable operators [7]. With this aim, we use the usual 
creation and annihilation operators respectively a} and a which are defined as follows 

+ p + imojx p — imujx 
a = /o 1 ' a = fo * ' (3) 

where p = —i-^- The external magnetic field is chosen in z-direction (i.e B = Bqz) 
in order to reduce the Hamiltonian defined in Eq.(l) and it has the form 

H = ^-a z + hucJa + p(a + a — cr_a^), (4) 

where e = 2pB . As a± = a T , it is pointed out that this Hamiltonian is not hermitian 

= ^a z + kioa^a — p(a + a — cr_a^), 

^ H. (5) 
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Thus as, 



PTH(PT)- 1 = -^a z + huja)a + p(a + a ] - <r_a), 



^ H, (6) 

one can see that the Hamiltonian (1) is not PT symmetric i.e H ^ H PT [1]. 

The next step is to write H in terms of differential operators(i.e p = —i4~) and 
of variable x. The purpose of these transformations is to reveal the exact solvability 
of the operator H by using the quasi-exactly solvable (QES) technique as has been 
considered in Ref. [14]. Replacing the operators at and a by their expressions(as given 
in Eq.(3))in the Eq.(4), the Hamiltonian of the model is written now as follows 

e p 2 — moj + m 2 u> 2 x 2 [a + (p — imoux) — a-(p + imujx)] . . 

M = a z+ 7T- + P rz = ("J 

In order to reveal the solvability of the above operator H, we first perform the stan- 
dard (often called " gauge") transformation 

H = R~ 1 HR , R = exp(-^-). (8) 

After some algebra, the new Hamiltonian H is obtained and is given by 

~ e Id 2 d W+p — O-iv + 2imuix)} 
H =-a z — 7 + ux— + p v —^- v 

2 2m dx 2 dx \ / 2mcuh 

e p 2 d [<r+p- <r-(p+ 2imux)} 

= —a z H V ojx- — hp- ; 9 

2 2m dx H V2m^Jh V ' 

Replacing a z , a + and cr_ by their matrix form, the final form of the Hamiltonian H 
reads 



H 



\ . 0-+ujx4z \ / 



2m^~ UJX dx U _ J y/2^JE 



(10) 



2™. +^^ + 2 ^ V / 2W[ 



Then, the operator /J is typically QES because it preserves a finite dimensional vector 
spaces of polynomials namely V n = (P n _i(x), P n (x)) 1 with n G N. Moreover # is 
exactly solvable because n does not have to be fixed (it can be any nonnegative 
integer) . 

Note that the above Hamiltonian H is not invariant under simultaneous parity 
operator(P) and time reversal (T)reflection (i.e respectively x — > —x and i — > —i) 
[1]. Even if the operator i/(therefore if)is nonhermitian and not PT invariant, it 
was pointed out that its spectrum is real. The reality of the eigenvalues of H is a 
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consequence of the unbroken Pa z (i.e combined parity operator P and Pauli matrice 
a z )invariance of H (i.e [H,Pa z ] = 0). In other words, the spectrum is real because 
H is pseudo-hermitian with respect to <r 2 (i.e o z Pl<j~ x = W) and also to the parity 
operator P (i.e PHP^ 1 = H^) [4, 5, 6]. We would like to mention that it is not 
necessary to calculate the energy eigenvalues and their corresponding eigenvectors 
of H because they have been determined in [5]. In the following section, we will 
construct the spectrum of the generalized Hamiltonian of the one given by Eq.(l). 



3 Family of exactly solvable Hamiltonians 

The original JCM is defined by the Hamiltonian 

H = -<r 3 + huaJa + p(acr + + aV_), (11) 

where p is a real parameter (i.e it is a real coupling constant). Note here that the 
Hamiltonian H is hermitian. 

In the next, we consider an extension of the above JCM Hamiltonian in the form 

H = ^o- 3 + Hua'a + P(a f a) + p(a k a + + 0(a t )V_), (12) 

where = ±1 and P(a^a) denotes a polynomial of degree d > 2, k is an integer > 1 
and p is an arbitrary real parameter. In fact, the above Hamiltonian is nonhermitian( 
i.e for = — 1 ) and not PT invariant but it satisfies the pseudo-hermiticity with the 
operators P (operator of parity) and 03 (Pauli matrice ). considering = +1 the 
Hamiltonian given by the Eq.(12) becomes hermitian. Both for these cases, it can 
be easily observed that the energy spectrum is entirely real. Thus, notice that the 
above Hamiltonian (12) is a generalization of the Hamiltonians given by the Eqs.(l) 
and (11). The matrix form of H reads 



/ %uoa)a + P(a)a) + | pa k 

V (f)p{a)) k hua)a + P{a)a) - 

which can be easily checked to preserve the vector spaces 



(13) 



V n = span{ ( "M.f „° , ■. )) , // ( IX. (1 !) 



J'\\n + k) 

It means that the action of the operator H on the vectors states ^ g ^ an d ^ | + fc\ ) 

can expressed as linear combinations of these same states. Here, we are allowed to 
conclude that H is exactly solvable because it preserves the vector space V„ for any 
integer n. 

The next step is to find the energy eigenvalues and their corresponding eigenvec- 
tors of the Hamiltonian H for = —1 and for = +1. For this purpose we recall the 
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following identities [5] 



a)a | n, \m s ) — n \ n, ^m s ), 

I 1 \ I 1 \ 

o- 3 I n, im s ) =m a \n, -m s ), 

IN lx lx 

0+ I n, 5 ) = ; ff + I n '^2^ = ' n ' 2'' 

1 1 

<7_|ra,-i) =0 ; cr_ I n,-) =\ n,--), (15) 

with n and m s = ±1 are respectively the eigenvalues of the number operator a^a 
and the operator a 3 . It is readily seen that the state | 0, — |) is a ground state 
of the operator if (i.e it is constructed by the lowest values of n and m s which are 
respectively and —1). We have now to consider the action of H to the state | 0, — |) 
in order to find its associated eigenvalue 

H I 0, -\) = -(7 3 I 0, --) + hwcta I 0, --) + P(ata) | 0, --) 
+pa k a+ I 0, ~) + 0p(at)V_ | 0, ~), 

= 2 <Ts| °'"2 > ' 

= -|l°'-|>- (16) 

It is proved now that — | is the eigenvalue of the ground state | 0, — |). It is easily 
understood that the next state | 0, |) is not an eigenstate alone of the Hamiltonian 
if because applying this operator to this state, we obtain a linear combination of two 
states I 0, |) and | k, — |), 

^|0,i> = ||0,l>±pVfc!|fc,~>. (17) 

The state | A;, — |) under the action of H leads to a linear combination also of two 
above states 

H I k, -\) = (huk + P{k) - E -) I k, -\) + P Vk\ I 0, 1>. (18) 

The excited states | 0, |) and | A;, — |) span an invariant subspace in the space of 
states so that the Hamiltonian matrix is written as follows 

/ e 



Hk= [(t ) pV¥. huk + P{k)-i) (19) 

In particular, note that for k — 1, P(fc) = 0(i.e P(fc) = A; d , d > 2) and considering 
= — 1, iffc becomes the matrix ifi constructed in [5]. In order to find the eigenvalues 
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of the Hamiltonian matrix(19), we have to solve the following usual equation(i.e 
characteristic polynomial) 



det(H k - All) = 0, 

(f-A P Vk\ \ 

\<ppVk\ huk + P{k)-^-X) 
AX 2 -4(huk + P(k))X + 2(huk + P(k))t-e 2 + <f)4k\p 2 =0. (20) 

After some algebra, the energy eigenvalues(i.e square-roots of the above equation in 
A) of H k are 



huk + P(k) + ^(hujk + P{k) - e) 2 + <f>4k\p 2 



k 



u huk + P(k) - ^(huk + P(k) - e) 2 + (j)4k\p 2 
— 2 ' 

It is easily checked that for k — 1, P(k) = and for <fi = —1), we obtain the eigenvalues 
A!' 77 determined in[5]. These are the energy eigenvalues of the Hamiltonian (1). The 
next step now is to calculate the associated eigenvectors of the above eigenvalues X{ ,n ■ 
Here, we propose to consider two cases : the first case for = —1 and the second one 
for = +1. 

3.1 The case = -1 

Considering <fi = —1, the eigenvalues (21) are given by 



I Tiuok + P(k) + yj (huok + P(k) - e) 2 - 4k\p 2 
h = 2 

huk + P(k) - J(hujk + P(k)-e) 2 -4k\p 2 

K 1 = v 2 — — ■ (22) 

For the sake simplicity, we can impose P{k) = and the eigenvalues X 1 / 1 have the 
form 



j hook + yj(hujk - e) 2 - 4k\p 2 
Xi, 



k 2 



huk - J (huk - e) 2 - 4k\p 2 
K 1 = Vl 2 " (23) 



The following relations are considered as in [5] 



\huk — e\ > 2p\f~k\, 
2pVk\ =(huk-e)sm9 k (24) 
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and the Hamiltonian matrix given by (19) reads 



e 



— pv^I frw/c — | / ' 

| \{hook — e) sin 

— \{%uik — e) sin /kuA; — | 

Taking account of the following equation 



(25) 



2 2 * 



!(fiwA;-e)sin0 fc \ M\ _ ,i,ii (A 



K' 11 ( d ) ' (26) 



-±(ftu;£;-e)sin0 fe hwk - % J \B J k \B 
the associated eigenvectors of X 1 / 1 are determined 

I Vk) = siny | 0,-) +cosy | k,--), 

for A£ = ^(l + cos 9 k )- | cos fc , (27) 

with A = sin ^ and £? = cos ^. 

I -0fe Z > = c os y | 0, i) + sin y | fc, -i), 

for Xi 1 = 5^(1 - cos0 fc ) + | cos0 fe , (28) 

with A = cos ^ and 5 = sin ^. 

In particular, for k — 1, it is easily checked that ipl and ipjf become respectively 
ip[ and ipl 1 which were determined in [5]. 

3.2 The case = +1 

Taking account of = +1 and imposing P(k) = , the eigenvalues (21) read 

l hook + ^{hook - e) 2 + 4k\p 2 

K = n > 



2 

The relations considered in Eq.(24) become 



hook - J (hook - e) 2 + Ak\p 2 
K 1 = V -■ (29) 



\hook-e\ > 2pVW, 
2pVk\ = {hook - e) sinh# fe . (30) 
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Following the same method used in the previous case, the eigenvectors associated to 
above eigenvalues (29) are written as follows 

l^fc) = sinhy | 0, i) + coshy | k,~), 

for \{ = ^^(1 + cosh 6 k ) - | cosh 9 k , 

I ibU) = cosh I 0, -) — sinh I k, — -), 
2 2 2 2 

/or A jr/ = ^(1 - cosh 4) + | cosh^ fc , (31) 

For if ^ if''' (i.e for = —1), it may be easily observed that two states given in 
(27) and (28) are not orthogonal to each other. But one can prove that the states given 
by Eq.(31) (i.e for = +1, H = W) are orthogonal. This property is a consequence of 
the hermiticity of H. In order to find the next excited states, one has to consider the 
next invariant subspace which is spanned by the vectors | 1, |) and | k + 1, — |). The 
eigenvalues and eigenvectors for this doublet can be determined following the same 
method used previously. 



3.3 The excited states 

The next step is to generalize the previous results to the invariant subspace which is 
spanned by the vectors | n, |) and | n + k, — |). Following the same technique used 
in the previous section and after some algebra, the Hamiltonian matrix for the above 
doublet is written as, 



/ hmn + P(n) + f p\Jn + l...y/n + k \ , . 

n+k ~ \<l>py/K+l...y/n + k hu(n + k) + P(n + k)-^J [ ' 

For the sake simplicity, we impose P(n) = P(n + k) = and H n+k is of the form 

tt _( hujn +^ pVn+ l...y/n + k\ , . 

n+k \<j>pyfr+l...yfc+k hw(n + k)-l J 1 j 

and its eigenvalues are 

j hu{2n + k) + ^Jihujk - e) 2 + 04p 2 (n + 1) ...(n + k) 



hu(2n + k) - ^/(htuk - e) 2 + 04p 2 (n + 1) . . . (n + k) 



KU = — ~ 2 — (34) 

In particular, putting k = 1 and = — 1 only in (34), the above eigenvalues be- 
come the eigenvalues A^+{ associated to the operator H given by the Eq.(l). These 
eigenvalues were determined in [5]. 
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Now putting 2py/n + 1 . . . \Jn + k = (htuk — e) sm9 n+k and 2py/n + 1 . . . Vn + k = 
(Hook — e) sinh.9 n+k in Eq.(34) respectively for = — 1 and for = +1, we find the 
eigenvectors corresponding to the doublet | n, |) and \n + k, — \). 

First considering = — 1, the eigenvectors associated to this doublet are 



■ hcos — 



V n+ k) = sin — | n, -> + cos^- | n + k,--), 



2 ' '2' ' 2 1 ' 2' 
/or A^ t+fe = ftum + -77- (1 + cos6> n+fc ) - ^ cos9 n+k , 



2 v ' 2 

I 1 II \ &n+k I 1 , . #ri+fc , , 1 > 

/or A^ +fe = Tiun + ^^(1 - cos n+k ) + | cos 6» n+fc . (35) 

Finally considering = +1 for the Eq.(34), the eigenvectors for the doublet | n, \) 
and | n + k, — \) are of the form 

| i^i+k) = sinh ^ | n, i) + cosh ^ | n + fc, ~), 

for A^ +fc = ftwn + ^7^(1 + cosh6» n+fe ) - | cosh6» n+A; , 

| V'n+ifc) = c °sh ^ | n, i) - sinh ^ | n + A;, ~), 

for A^ +fc = hun + ^^(1 - cosh6» n+fc ) + | cosh0 n+fe . (36) 

Note that all the discussions considered in the previous section are confirmed by these 
generalized results. 

4 Quasi-exactly solvable Hamiltonians 

In this section let us consider an extension of the Jaynes-Cummings Hamiltonian 
which includes two-photon interaction 

H 2 = -03 + hua^a + p(a + a 2 + a^a) 2 ) (37) 

The matrix form of the above Hamiltonian reads 

„ / huata + % pa 2 \ , . 

^ 2 = { p(at) 2 Wa - f ) ■ (38) 

It is clear that this Hamiltonian H is similar of the one reported in Ref.[ll] and is 
also a particular case of the Hamiltonian given in Eq.(13) (i.e if k = 2, P(a^a) = 0) 
and one can prove easily its exact solvability. However, if one would like to construct 
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an JC-type Hamiltonian including both a one-photon and a two-photon interaction, 
the above Hamiltonian should be modified as follows 



H 



12 



Tiuja)a + 



| per + pxa 

€_ 

2 



(39) 



where p, p 1 , pi are, a priori, arbitrary constants. 

Unfortunately, the corresponding operator if 12 is not anylonger exactly solvable. 
Indeed, it is easy to show that it fails to admit any finite dimensional invariant vector 
spaces. Accordingly, it is impossible (to our knowledge) to find its energy spectrum 
by algebraic methods. 

In order to restaure, at least partly, a certain algebraic solvability of Hi 2 , one 
can attemp to supplement the Hamitonian H\ 2 with an appropriate interation term. 
After some algebra, one can convince oneself that adding an interaction term of the 
form 



H, = 





n \p\(Jaat 



piacJa 




(40) 



H 7 



(41) 



leads to a new Hamiltonian Ht = H\ 2 + Hi which is quasi-exactly solvable, as we 
will now demonstrate. 

Assuming n to be an integer and redefining c = — c = — the operator H T 
reads 

huja'a + | pa 2 + ca(a'a — n) 

, 4>p{a)) 2 + c(a^a — n)a) huaJa — | 

where that a* and a are respectively the usual creation and annihilation operator and 
e is chosen as previously according to e = 2pBo- 

The main idea now is to reveal that the above operator Ht is quasi-exactly solv- 
able(QES). In this purpose we construct a finite dimensional vector space which is 
invariant under the action of H T - Let us apply now the Hamiltonian H to the states 



JV) 




and 





M) 



with JV.IeNas follows 

(RuJV + §) I JV) + Py J M(M - 1) | M - 2) + cVM(M - n) | M - 1) 
<t>p^/(N + 1)(JV + 2) | JV + 2) + c^/WTl(N + 1 - n) | JV + 1) + (RwM - f ) | M) 



(42) 



In order to be in agreement with the invariance of the two vectors states 




JV) 




and 



M) 



under the action of the Hamiltonian Ht, we have to impose the value 

of the integer n according ton = M = N + 2 (i.e N = M — 2). Taking account of 
the above fixed value of n, we obtain 



(KwN + f ) + p^/(JV + 2)(JV + l) 



I N) 



V 



I JV + 2) - cVJV + 1 | JV + 1) 



(43) 



ft^(JV + 2) - § + 4>P^(N+ l)(JV + 2) 

Finally the Hamiltonian ify is of the new form 

fiWa+§ pa 2 + ca(a^a- (JV + 2)) 



#7 



±p(at)2 + c( a ta — (JV + 2))at 
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(44) 



As it is clear from the Eq.(43), the Hamiltonian preserves the finite dimensional 
vector space namely 

V^spanjC^),^) , j = N,...,0; fc = iV + 2, . . . , o|, (45) 

and n is fixed according to n = N + 2. From this, we conclude that the Hamiltonian 
H T is quasi-exactly solvable. Hence the terms of perturbation added to H 12 have 
broken its non solvability. 

Notice that is also easily to reveal the quasi-exact solvability of the operator 
expressed in Eq.(41) by considering the matrix Hamiltonian Eq.(41) in terms of dif- 
ferential expressions. Here we have to replace the operators a' and a respectively by 
their differential expressions given by Eq.(3), performing the standard gauge trans- 
formation as, 

H T = exp(^-) H T exp(-^-), (46) 

and thus, after some algebra, we obtain a matrix Hamiltonian which preserves the 
finite dimensional vector space of the form Vk = (Pk(x), Pk+2{x)) t with k G N and 
n = k + 2 (i.en which is expressed in Eq.(41)). This operator H T (therefore H T ) is 
quasi-exactly solvable because it is expressed in terms of the integer n which is fixed 
according to n = k + 2. 



5 Spectral properties 

In this section, we would like to emphasize a few properties of the spectrum of the 
Hamiltonian discussed above. First we stress that for given k the JC model admits 
k levels which are p-independant and which are not involved in the list given above. 
They are of the form 




< j < k - 1 , 



where denotes the null vector of the Hilbert space. The corresponding eigenvalue 
is Ej = j - f . 

The spectrum of the JC model (and of its generalisations for k > 1) varies con- 
siderably with the parameter p. In Fig. 1, we show the evolution of six levels in the 
k — 2, (f) — 1 case. They correspond to the two p-independant eigenstates and the 
ones with n — 0, 1 in Eq.(34). In Fig. 1 and in the following we assume e = 1 for 
simplicity but the features pointed out below remain similar for e ^ 1. The same 
levels corresponding to the non hermitian case = — 1 are reported on Fig. 2. The 
contrast with Fig.l is obvious. Couples of eigenvalues regularly disappear at finite 
values of the coupling constants p. So that, at finite p only a finite number of real 
eigenvalues subsist, the other being real. In this respect, the Hamiltonian is like a 
quasi exactly solvable operator. 

The energy levels displayed on Fig.l corresponds to the six lowest ones in the 
limit p = 0. The figure clearly shows that they mix relatively quickly for increasing p 
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and that, for instance, eigenvectors involving two or more quanta become the ground 
state for p ~ 1. 

We have studied the evolution of the spectrum when the QES-extension of the 
model, if 12 = pa 2 +9a{l — j^a^a) namely characterized by the new coupling constant 

9, is progressivel switched on. Notice that the vector ip = (0 , |0))' is an eigenvector 
with E = — e/2, irrespectively of p, 9 

In the case p — 0, N — 1 the effect of the new term on the eigenvalues under 
consideration leads to 

E = ~\, ^(3 ±40) , i(9±2v/2>, \ 

These levels are indicated on Fig. 3 by the dotted lines and it is clearly seen that 
they also lead to numerous level crossing. 

The evolution of the eigenvalues corresponding to the case p = 1 is displayed 
by the dashed lines in Fig.3, supplemented by the black line E — —1/2 which is 
present irrespectively of p. The figure clearly shows that the occurence of the new 
term induced only one level mixing, namely two levels cross at E — —1/2 for 9 = 1.5 
For larger values of p, e.g. p — 2, the analysis reveals that the algebraic eigenvalues 
depend only weakly of 9. 



6 Series expansion and Recurence relations 

Here we would like to present another aspect of the QES Hamiltonian presented in 
the previous section. Following the ideas of [15] we will construct the solution for 
energy E under the form of a formal serie in the basic vector whose coefficients are 
polynomials in E. More precisely, we write the solution of the equation 

H T iP = Eijj, (47) 

in the form 

i>=(^*.'£WL) (48) 



T.T--ME) I J +2) 

and where Ht is given by the Eq.(41). After some algebra it can be seen that the 
polynomials pj(E),qj(E) obey the following recurence relations 

A j+l P j+1 + BjPj = Q, (49) 

where 



Aj+i = ( P}J(j + 2)(j + 3) —(E - (j + 1)_ 



c(j + 2-n)VJT2, 

c(i + 2-n)VJT2 



(E-(j + 2) + l) pJ(j + l)(j + 2) 



Pj ={p ) > J = -2,-1,0,1,... (50) 
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These equations have to be solved with the initial conditions 



q- 2 = 0,q-i=M (51) 

with M fixing the normalisation of the solution. Then the solution for qj turns out 
to be a polynomial of degree E 2 K The quasi-exact solvability of the system leads to 
the fact that i„_i is not invertible and that p n -\ can be choosen arbitrarily. With 
the choice p n -\ = it turns out that all polynomials Pj,qj with j > n — 2 are 
proportional to q n ^(E). As a consequence for fixed n and for the values of E such 
that q n s(E) = the serie above is truncated and the set of algebraic eigenvectors are 
recovered. We would like to stress that series considered in this section are built with 
the basis vector of the harmonic oscillator and not on monomials in x contrasting 
with the construction of Ref.[15]. In the case of standard QES equations [15] there it 
appears a three terms recurence relations which leads to sets of orthogonal relation. 
In the case of systems of QES equations adressed in [16] the recurence relation is also 
three terms but the situation here is quite different. Actually it is to our knowledge, 
an open question to know whether the set of polynomials (pj(E), qj(E)) are somehow 
orthogonal as it is the case for standard scalar equation. 

7 Hidden algebraic structures 

As pointed out in the previous sections, the different Hamiltonians studied here posses 
the property that their spectrum can be (partly or fully) computed. This property 
is deeply related to the fact that the corresponding operators are elements of the 
enveloping algebra of particular graded algebra in an appropriate finite dimensional 
representation. The classification of linear operators preserving the vector spaces 
V(m, n) = (P m (x) , P n (x)Y was reported in [10]. It is shown that these operators are 
the elements of the enveloping algebra of some non-linear graded algebra depending 
essentially of \m — n\. Note that, in the present context, the difference \m — n\ is 
nothing else but the parameter called k in Sect. 3. The cases k — 1 and k — 2 are 
special because the underlying algebra is indeed a graded Lie algebra. In the case 
k — 1, related to the conventional JC model, the Hamiltonian is an element of the 
enveloping algebra of osp(2, 2); in the representation constucted in [9]. The generators 
involved in this relation do not depend explicitely on n, i.e. of the dimension of the 
representation, explaining that the Hamiltonian is exactly solvable. Finally, in the 
case k = 2, the Hamiltonian is an element of the graded Lie algebra q(2), as shown in 
[17, 18]. This algebra possesses an sl(2)xU(l) bosonic subalgebra and six fermionic 
operators splitted into three triplets of the sl(2) subalgebra. In the case of the JC 
model corresponding to k — 2, the Hamiltonian is independant on the dimension of 
the representation n and the model is exactly solvable. For the modified model of 
Sect. 4, the supplementary interaction term Hi defined in (40) indeed depends on 
n and the operator admit only the vector space V n as finite dimensional invariant 
vector space. 
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8 Conclusions 



In this letter, we have considered several extensions of JCM by adding to its origi- 
nal Hamiltonian the polynomial P(a)a) of degree d > 2 and an arbitrary sign, say 
= ±1, in the non-diagonal interaction term. In fact, considering this sign = — 1, 
these extended Hamiltonians are nonhermitian and not PT invariant but they sat- 
isfy the pseudo-hermit icity with respect of different operators P and a 3 . This new 
property reveals the reality of the energy spectrum which has been constructed alge- 
braically. They become hermitian when one considers the sign = 1. Notice that 
these Hamiltonians are completely solvable as it has been pointed out by the QES 
technique. 

Several usual properties available with hermitian Hamiltonian are not kep with 
pseudo-hermitican. Namely the eigenstates given by (27) and (28)(i.e corresponding 
to the doublet | 0, |) and | k, — |)) are not orthogonal to each other, but they are 
orthogonal to all eigenstates corresponding to other doublets. For example, the eigen- 
state (27) and the one given by Eq.(35)(i.e it corresponds to the doublet | n, |) and 
| n + k, — |)) are orthogonal to each other. The eigenstates of any particular doublet 
are orthogonal to each other only if 9 m = mix (i.e with m — 0, 1, 2, . . . , k, . . . , n + k), 
this implies p = because it depends to sin6* m . In fact, as the energy eigenvalues are 
entirely real, it is impossible to have all eigenstates orthogonal to each other. This is 
explained by the unbroken symmetry of the operator Pcr 3 . But for energy eigenvalues 
complex, the orthonormality condition is satisfied by all the associated eigenstates. 
All these discussions are the result of the scalar product applied to those eigenstates. 

We manage to construct a JC-type Hamiltonian describing both one and two- 
photons interactions in terms of quasi exactly solvable operators. This involves a 
very specific interaction term of degree one in the creators and annihilators which can 
be seen as a perturbation of more conventional p-photons interacting term. Several 
properties of this new family of QES-operators have been presented. Namely, (i) 
they can be written in terms of the generators of the graded Lie algebra osp(2,2) 
in a suitable representation; (ii) when expressed as series, the formal solutions of 
Ht^P = Eip leads to a different type of recurence relation between the different terms 
of the series. 
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Figure 1: The first few energy levels in the k = 2-JC Hamiltonian for e = 1 and 
= 1. 
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Figure 2: The first few energy levels in the k = 2-JC Hamiltonian for e = 1 and 
6 = -1. 
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Figure 3: The first few energy levels in the QES deformed k = 2 JC Hamiltonian as 
function of the parameter 9, the energy level E = —1/2 (solid line) is independant of 
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